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As one of the brightest active blazars in both X-ray and very high energy 7-ray 
bands, Mrk 501 is very useful for physics associated with jets from AGNs. The 
ARGO-YBJ experiment is monitoring it for 7-rays above 0.3 TeV since November 
2007. Starting from October 2011 the largest flare since 2005 is observed, which 
lasts to about April 2012. In this paper, a detailed analysis is reported. During 
the brightest 7-rays flaring episodes from October 17 to November 22, 2011, an 
excess of the event rate over 6 a is detected by ARGO-YBJ in the direction 
of Mrk 501, corresponding to an increase of the 7-ray flux above 1 TeV by a 
factor of 6.6±2.2 from its steady emission. In particular, the 7-ray flux above 8 
TeV is detected with a significance better than 4 a. Based on time-dependent 
synchrotron self-Compton (SSC) processes, the broad-band energy spectrum is 
interpreted as the emission from an electron energy distribution parameterized 
with a single power-law function with an exponential cutoff at its high energy 
end. The average spectral energy distribution for the steady emission is well 
described by this simple one-zone SSC model. However, the detection of 7-rays 
above 8 TeV during the flare challenges this model due to the hardness of the 
spectra. Correlations between X-rays and 7-rays are also investigated. 

Subject headings: gamma rays: general - BL Lacertae objects: individual (Markar- 
ian 501) - galaxies: active - radiation mechanisms: non-thermal 



Introduction 



Blazars, including BL Lac objects and flat-spectrum radio quasars (FSRQs), are the 
most extreme subclass of active galactic nuclei (AGN). Most of the identified extragalac- 
tic 7-ray sources belong to this category. Their emission is believed to be dominated by 
non-thermal and strongly Doppler-boosted radiation from a relativistic jet of magnetized 
plasma which is aligned along our line of sight. The physical mechanism for the production 
of their 7-ray emission is still under debate. The leptonic models attribute the 7-ray emission 
to the inverse Compton scattering of the synchrotron (synchrotron self-Compton, SSC) or 
external photons (external Compton, EC ) by the same pop ulation of relativistic electrons 
(jGhisellini et al.lll998t iDermer et al.lll992l ; ISikora et al.lll994f ). therefore an X-ray/7-ray cor- 
relation is expected. The lack of strong emission lines in the radiation from BL Lac objects is 



taken as one of the evidences for a minor role of ambient photons (e.g.. iKrawczynski I (12004 )). 
and hence the SSC model is favored. The hadronic models attribute the 7-ray emission to 
prot on-initiated casca des and/or proton-synchrotron emission in a magnetic field-dominated 



jet (lAharonian 1 12000| ). However, the tight X-ray and very high energy (VHE) 7-ray cor- 
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relation and the very rapid 7-ray variability are taken as the strong challenges to models 
based on hadronic processes. Recently, a long-term continuous monitoring of Mrk 421 has 
been performed base d on the ARGO-YBJ experiment and satellite borne X-ray detectors 
(jBartoli et al.l 1201 laj ) . According to this investigation, both the temporal and the spectral 
results generally favor the SSC model. Even in the framework of the SSC model, the fun- 
damental question referred to the origin of the flux and spectral variability, observed on 
timescales from minutes to tens of years, is still open. 

Mrk 501 (z = 0.034) was discovered with VHE emission by the Whipple collaboration 
(jQuinn et al.lll996l ). It is one of the best-studied blazars with extensive studies on various 



timescales. In 1997, Mrk 501 went into a state with surprisingly high activity and strong 
vari ability and became more than a factor 10 brighter (above 1 TeV) than t he Crab Neb- 



ula (IDjannati-Atai et al.l Il999i ; lAharonian et al.l Il999l ; lAmenomori et al 



7-ray flux variability on a timescale of minutes was observed in 2005 ({Albert et al.l 120071 ). 



•200(3). The fastest 



Significant spectral variability was d etected with the harder spectrum at the brighter states 



comp ared to the low- activity states (lAlbert et al.l 120071 ; lAnderhub et al.ll2009l ; lAcciari et al. 



201ll ). In the middle of 2009, a multi-frequency observational campaign for 4.5 months on 
Mrk 501 was carried o ut with excellent e nergy coverage from radio to VHE 7-ray when it un- 
derwent a low activity (lAbdo et aU201ll ). Throughout the campaign, the source was sampled 
quite uniformly in all wavelength bands except for the VHE band, because the Cherenkov 
Telescopes cannot operate during non-optimal weather conditions or bright moonlight peri- 
ods. In October 2011, Mrk 501 underwent a strong flare detected by the MA XI satellite in 



X-rays (jSootome et al.ll201ll ) and by the ARGO-YBJ detector in VHE 7-rays (jBartoli et al. 
20 lid ). 



To understand the variability of emission and the underlying acceleration and radiation 
mechanisms in jets, continuous multi-wavelength observations from the X-ray to the VHE 
7-ray band are crucial especially over a very long term. The broad band energy spectra 
could provide constraints on the parameters of the models. The Cherenkov telescopes can- 
not constantly monitor AGNs because of their limited duty cycle and narrow field of view 
(FOV). The wide-FOV ARGO-YBJ detector, operated with high duty cycle (> 85%), is 
more suitable for monitoring. Working at energies above 300 GeV, ARGO-YBJ extends 
the multi-wavelength survey carried out by the satellite borne X-ray detector Swift and the 
GeV 7-ray detector Fermi-LAT. Particularly, the spectral energy distribution (SED) of Mrk 
501 is covered without any gap from 100 MeV to 10 TeV. All the measurements would set 
strong constrains on the model of the emission of AGNs. 



In this paper, we report on the multi- wavelength view of the emission from Mrk 501 from 
August 2008 to April 2012, including the average spectra during quasi-steady and flaring 
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periods. 



2. The ARGO-YBJ experiment and data analysis 



The ARGO-YBJ experiment, located at the Cosmic Ray Laboratory of Yangbajing 
(Tibet, P.R. China) at an altitude of 4300 m a.s.L, is the result of a collaboration among 
Chinese and Italian institutions and is designed for VHE 7-ray astronomy and cosmic ray 
observations. The detector consists of a single layer of Resistive Plate Chambers (RPCs), 
which are equipped with charge readout strips (6.75 cm x 61.80 cm each). The logical OR 
of signals from eight neighboring strips constitutes the pixel (called a "pad") for triggering 
and timing purposes. Each RPC is read with ten pixels. 130 clusters (each composed of 12 
RPCs) are installed to form a carpet of about 5600 m 2 with an active area of ~93%. This 
central carpet is surrounded by 23 additional clusters ("guard ring"). The total area of the 
array is 110 m x 100 m. The ARGOYBJ detector is operated by requiring the number of 
fired pads (N pai j) to be at least 20 within 420 ns on the entire carpet. The trigger rate is 
3.5 kHz with a dead time of 4%. The angular resolution, pointing accuracy and stability 
of the ARGO-YBJ detect or array have been t horoughly tested by measuring the shadow of 
the Moon in cosmic rays ( IBartoli et al.ll2011b| ). The absolute energy scale uncertainty is less 

More details about the 

313). 



than 13% for all measured cosmic ray showers (jBarto 



detector and the RPC performance can be found in (lAielli et al. 



i et al 



201 



b) 



2006, 2009a 



For the multi-wavelength investigation together with Fermi-LAT, the data collected 
by ARGO-YBJ after August 2008 when Fermi was launched are used. The total effective 
observation time is 1179.6 days. To achieve a good angular resolution, events with ze nith 
angle less than 50° are used, and further selection criteria ( IBartoli et al.l 1201 lal . |2012| ) are 
applied. The total number of events after filtering is 1.86 xlO 11 for this work. No 7/hadron 
discrimination is applied. The opening angle ^70, which contains 71.5% of the events from 
a point-like source, is 1.36° for events with N pac i > 60. In order to remove the effect of 



the cosmic-ray anisotropy, the method described in (IBartoli et al. 



2011ai 12012( 1 was applied. 



The significance of the excess is estimated using the Li fc Ma ( 19831 ) method. With this 
data analysis, the significance of the excess observed from the direction of the Crab Nebula 
is 17 standard deviations (er) in 3.5 years, which indicates that the 3.5 -year cumulative 5 a 



sensitivity of ARGO-YBJ has reached 0.3 Crab unit for point sources ( jCao fc Chen 1 1201 lh . 
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3. Results 
3.1. Light curves 

The daily flux from Mrk 501 at energy 15-50keV provided by Swift/BAT^ is publicly 
available and used in this work. The light curve from August 2008 to April 2012 is shown 
in panel (a) of Figure [1] with a bin size of 30 days. The best fit with a constant value 
for the light curve is (8.9±0.4) xlO -4 counts cm -2 s _1 with a x 2 °f 492.9 for 44 degrees 
of freedom (ndf). A significant feature is the flare at the end of 2011 with the flux being 
enhanced by a factor about 4. Without data during the flaring period, the best fit with 
a constant value for the light curve is (6.0±0.4) xlO -4 counts cm -2 s _1 with a x 2 / n df of 
69.3/38. This result could be an evidence for a small X-ray variability before this flaring 
period. Zooming in the Swift/BAT light curve at 15—50 keV shown in Figure |2J the large 
flare started in October 17, 2011 (MJD=55851) and decreased to a low-activity state around 
the average level on November 22 (MJD=55887) (flare 1 hereafter). Afterwards, Mrk 501 
became increasingly active for a longer period up to about April 2012. Its brightest flaring 
episode is in November 8, 2011, during the flare 1 period. 



The Fermi-LAT data were analyzed using the ScienceToolsj. The light curve is gener- 
ated using aperture photometry. The panel (b) of Figure [1] shows the flux at energies greater 
than 0.3 GeV that is contained within a 2-degree cone centered on Mrk 501, which is the 68% 
containment angle of the reconstructed incoming photon direction for normal incidence. A fit 
with a constant value yields a x 2 / n df of 107.7/44 indicating a moderately variable behavior 
consistent wit h the X-ray analys i s. Th e Discrete Correlation Function (DCF), computed as 



prescribed by lEdelson fc Krolikl ( 119881 ). for the BAT/LAT data points shown in Figured] is 



DCF = 0.63± 0.26 for a time lag of zero, which is greater than the previously measured DCF 



= 0.32±0.22 ( lAbdo et al.ll201lf ). Since the significance in both analysis is less than 2.5cx, 
only minor correlations are observed. During the X-ray flaring period, the GeV 7-ray flux 
increased above the long-term average, but this flux increase is not significant. Therefore, 
the light curve does not indicate a significant correlation with the X-ray data during the 
flare. 

The light curve in the TeV 7-ray range detected by ARGO-YBJ is shown in panel (c) of 
Figure [TJ During the X-ray flare, the flux of TeV 7-rays also increases. A fit with a constant 
emitting rate yields a x 2 / n df of 71.9/44, while the x 2 / n df is reduced to 42.59/38 by simply 



1 Transient monitor results provided by the Swift/BAT team: http://heasarc.gsfc.nasa.gov/docs/swift/ 
results / transients /weak/Mrk501/. 

2 http : / /fermi . gsfc . nasa.gov/ssc / 
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excluding the data during the X-ray flares. The TeV 7-ray flare 1 was clearly detected as a 
counterpart of that in the X-ray band, and the count rate increase is by a factor about 6. 
The DCF for the BAT /ARGO-YB J data points shown in Figured] is 0.85±0.36 for a time 
lag of zero, while the DCF for the L AT / ARGO-YB J is 0.44±0.31. 



3.2. Photon Energy Spectra 

To investigate the evolution of the spectra, their time-averages during the long-term 
quasi-steady state from August 5, 2008 to October 16, 2011 and during flare 1 are estimated 
separately. The integrated flux from Mrk 501 observed by ARGO-YBJ has a statistical 
significance of 5 o (see panel (a) of Figure [3]). During flare 1, the flux from Mrk 501 was 
detected by ARGO-YBJ with 6.1 a significance (see panel (b) of Figure |3]), corresponding to 
an increase of the 7-ray flux above 1 TeV by a factor of 6.6 from its steady emission. Events 
with Npad > 60 are used for both panels of Figure [3J For comparison, all the spectra are 



show n in FigureEJ where an average spectrum over 4.5 months was adapted from (lAbdo et al. 



20111 ). Without significant activity observed in the last 3 years before flare 1, the average 
spectra over 4.5 months could approximately represent the average spectra over 3 years. All 
the other spectra presented in Figure [6] are estimated as described in the following. 



3.2.1. Swift: X-ray 



The Swift /XRT is a focusing X-ray telescope with energy range from 0.2 to 10 keV. 
In Windowed Timing mode during flare 1, XRT data in four time windows are available 
from HEASARC |, and the exposure time is about 1 ks in each window. The XRT data set 
was processed with the XRTDAS software package (v. 2. 6.0) following the standard recom- 
mendations. The XRT average spectrum in the 0.5—10 keV energy band was fitted using 
the XSPEC package (v. 12. 7.0). We adopted a power-law model for the photon-flux spectral 
density, with an absorption hydrogen-equivalent column density fixed to the G alactic value 



in the direction of the source, namely 1.56 x 10 2U cm 2 ( iKalberla et al.l 120051 ). We obtain 



the spectrum of (0.0407 ± 0.0007) x (jg^) 



)20 

-1.824±0.022 



ph cm 



-2 g -l 



keV" 



3 http://heasarc. gsfc.nasa.gov/ 
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3.2.2. Fermi-LAT: HE 7-rays 

The LAT data from a region centered on Mrk 501 with radius of 10° was used to estimate 
the spectrum. The data analysis was performed following the standard recommendations by 
using an unbinned maximum-likelihood method to estimate the source SED. We adopted a 
power-law model for the energy range from 0.1 to 300 GeV. During the long-term period, 
Mrk 501 was detected with a Test Statistic (TS) value of 5118.2 (~ 71.5cr). We obtain the 
spectrum (2.300±0.068) x 10~ 12 ( n^iv^v) _1 - 7r4±0 021 phcm -2 s" 1 MeV -1 . This is consistent 
with the re sult reported by t he LAT collaboration using the data collected during the first 



24 months ( lAbdo et al.ll2012l ). The spectrum was found to be very stable in the GeV band 
in spite of a moderate variability over a long period. During flare 1, the TS was 329.1 (~ 
18.1a) and the spectrum was (3.53 ±0.46) x 10~ 12 ( 18 44.f Me y )~ 1,640±0 ' 084 ph cm" 2 s" 1 MeV" 1 . 
Compared with the long-term result, the flux increased slightly. 



3.2.3. ARGO-YBJ: VHE 1 -rays 
The VHE 7-ray spectrum was estimated using a distribution of the excess i n the number 



of even ts as a function of N pa d. A widely used procedure, that is described in (jBartoli et al. 



2011al ). was followed. In this procedure, the spectrum of Mrk 501 was assumed to be a 



power law. The ARGO-YBJ detector response has been taken into account using G4argo 



(IGuo et al.ll2010l ). The simulated events are sampled in the energy range from 10 GeV to 
100 TeV. We define six intervals with N pad of 20-59, 60-99, 100-199, 200-499, 500-999 
and >1000. The best fit gives a differential flux(cm~ 2 s _1 TeV -1 ) 

/ j? \ -2.59±0.27 

(1.92 ± 0.44) x 10- 12 ( — J (1) 

for the long-term period, corresponding to 0.312±0.076 Crab units above 1 TeV. The median 
energies of the six intervals are 0.45, 0.89, 1.4, 2.8, 5.6 and 11 TeV, respectively. Both the 
flux and the spectral index are similar t o those obtained b y VERITAS and MAGIC during 



the low-activity state, as presented in (lAbdo et al.l 1201 ll ). In particular, the ARGO-YBJ 



data show a smooth extension of the Fermi-LAT results, as evident in Figure EJ During 
flare 1, the differential flux (cm -2 s -1 TeV -1 ) is 



-2.07±0.21 

1-12 



(2.92 ± 0.52) xlO^f^^j (2) 



corresponding to 2.05±0.48 Crab units above 1 TeV, that is a factor 6.6±2.2 compared with 
its long-term steady state. The median energies of the six intervals are 0.89, 1.1, 1.8, 3.5, 



-9 - 



7.1 and 14 TeV, respectively. Only the statistical error is qu oted here, and the systematic 
uncertainty in the flux measurement is estimated to be <30% lAielli et al.l ( 120101 ) . 



Discussion 



4.1. Spectra Corrected for the EBL Absorption 

During flare 1, 7-rays with a median energy of 8.4 TeV are also observed with a signifi- 
cance greater than 4cr, as shown in Figure HJ 7-rays with such a high energy from Mrk 501 
have not been detected since the 1997 flare. The SEP at energies ab ove 0.9 TeV is harder 
than those observed dur ing the flares in 1997 (lAharonian et al.ll200ll ) and in June 30, 2005 
flare (lAlbert et al.ll2007l ). as shown in Figure El although the spectral indices are consistent if 
the statistical error is taken into account. Mrk 501 is a nearby source, so we do not expect a 
significant absorption of its intrinsic source spectrum due to Extragalactic Background Light 
(EBL) at energies below 1 TeV, while the absorption at higher energies is still considerable. 
Therefore, it is useful to test different EBL models assuming a minimum intrinsic photon 
spectral index. A natural minimum spectral index is 1.64, constrained by the spectrum in 
the GeV band, since that at higher energies should be steeper. 

Here we use four kinds of models with different flux levels of the EBL, among many 
models. Assuming a single power law for the VHE flux, the indices of the derived unab- 
sorbed spectra using different models are:(l)l. 801^29 f° r the "low-IR" model proposed by 



EBL model of 



Kneiske et al.l f 20041). which giv e s a si milar result co mpared with that o btained using the 



Aharonian et al.l (120061) . according to lAlbert et al.l (120071 ): (2) 1.45 



; t°o2 for 



the model of iFranceschini et al.l ( 120081 ). which is widely used to correct the VHE SED of 



extragalactic sources; (3,4) 1 -30^ and 1.11 

-oil f° r the baseline and fast evolution models 



proposed by IStecker et al.l (120061 . 120071 ). respectively. For comparison, all the unabsorbed 
spectra are shown in Figure EJ The spectral indices obtained using models (2), (3) and 
(4) exceed the minimum spectral index boundary of 1.64, however they are consistent with 
this limit if the statistical error is taken into account. It has to be noted that models (3) 
and (4) have been excluded w i th higher significance by pr evious tests carried out around 



1 TeV (lAharonian et al.l l2006t iGeorganopoulos et al.l 120101 ) and tens of GeV ( lAbdo et al. 



2009L boioh . Since our data extend to about 10 TeV, the corresponding EBL photon energy 
is substantially lower. The EBL model with the minimum absorption is used when modeling 
the SED in the following section. 
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4.2. Modeling of the overall SED 



The long-term averaged SED, especially the continuous measurement of the second 
component in the energy range 0.1 GeV-10 TeV obtained in this work, provides us a robust 
baseline for an insight into the underlying physics in Mrk 501. Over the very wide radio- 
VHE ene rgy range, we fitted a one-zone SSC model pro posed by iMastichiadis &: Kirkl (119951 ) 
(see also IMastichiadis fc Kirk! 119971 ; I Yang et al.l 120081 ) to the SED. There are a few free 
parameters to be determined in this model, including the Doppler factor 5 = l/[r(l — 
(3 cos 8)}, the spherical blob radius R, the magnetic field strength B, the electron spectral 
index s, the electron maximum Lorentz factor 7 max , and the electron injection compactness 
l e = ^rrieCa^R 2 J 1 °°rf7(7 — l)Q e - in order to determine the Doppler factor and the injection 
compactness, further parameters must be determined, i.e., the Lorentz factor T, the speed of 
the blob c/3 and the Lorentz factor of electrons 7. Moreover, o"t is the Thomson cross section, 
9 the angle between the direction of motion of the blob and the observer's line of sight, and 
Q e the electron spectrum at injection, which is assumed to be a power law cutting off at 
7max, with a normalization factor q e , Q e = q e ^~ s exp(— 7/7 max )- This model is different from 
other one-zone SSC models which introduce more fr ee parameters using stationary injected 



electr on spectra with a double power-l aw function (jTavecchio et al.l l200ll ; lAnderhub et al. 



20091 ) or even triple power-law function (lAbdo et al.ll201ll ) with an exponential cutoff at their 
high energy end. A full time-dependent evolution of the electron and photon spectra based 
on this mod el was simulated for a given injection spectrum of electrons. The low-IR model 
proposed by iKneiske et al.l ( 120041 ) is used to take into account the absorption of 7-rays in 
the EBL when modeling the SED. 

The best fit to the long-term SED is shown in Figure El with the corresponding param- 
eters given in Table 1. General agreement between the model and the data is achieved with 
this simple one-zone SSC model. The parameters in the model are found in general agree 



ment with those found in the previous analysis on the similar source Mrk 421 ( iBartoli et al. 



2011al ). The best fit for the flare 1 is also shown in Figure [6j with the corresponding param- 
eters given in Table 1. The highest-energy data points of the SED (> 6 TeV) during flare 
1 cannot be well reproduced by only modifying the parameters. It is important to point 
out that the X-ray spectrum becomes harder during the flare. The peak energy is shifted 
to about 10 keV during flare 1 from about 1 keV in the quasi-steady state. The shift of the 
X/7-ray peak to h igher energies during flares is a common f eature which has b een reported 
many times (e.g. (lAlbert et al.l 120071 ; lAnderhub et al.l 120091 ; lAcciari et al.l 1201 ih ) . However, 
the detection of >6 TeV 7-rays with an energy flux similar to that at 1 TeV is unusual. In the 
framework of the SSC model, it is difficult to reproduce such high energy 7-rays as detected 
by ARGO-YBJ. Since the 7-rays above 1 TeV are typically produced in the Klein-Nishina 
regime, their rate should be strongly suppressed. The radiation mechanism during flares 



may be different from that in the quasi-steady state. Differ ent radiation mechanisms , such 
as more complex SSC models or hadronic processes (e.g. in lAbramowski et al.l (120121 )). are 
needed to improve the understanding of the flaring phenomena. 



SUMMARY 



We have presented a continuous long-term monitoring of Mrk 501 from August 2008 
to April 2012. Both the flux and the spectral index are consistent with tho se obtained 



by V ERITAS and MAGIC during a 4.5 month long multifrequency campaign fjAbdo et al. 
20111 ). Combining the observations by the ARGO-YBJ experiment with the space-borne 



experiments Swift and Fermi, the investigation was performed over a wide ener gy range 
from 0.5 keV up to 10 TeV (a value higher than the maximum in (lAbdo et al.1120111 )). Using 
all the data covering various energy bands during the quasi-steady phase of the blazar, its 
SED is fitted with a simple one-zone SSC model assuming a single power-law with expo- 
nential cutoff for the electron spectrum at injection. The model parameters are found in 
agreement with those resulting from previous analyses for various AGNs, indicating that 
similar radiation mechanisms are in action. A strong flare of the blazar in the VHE region 
was observed by ARGO-YBJ in October 2011, while no Cherenkov Telescope could observe 
Mrk 501 during this period. It is well determined to be a counterpart of the X-ray flare in 
the same time period, giving for the BAT/ARGO-YBJ data points a DCF of 0.85±0.36 for 
a time lag of zero. On the contrary, there is no significant increase of the flux at energies 
around 1 GeV. Remarkably, 7-rays with energies above 8 TeV are detected, which did not 
happen since the 1997 flare. The spectral shape obtained wit h the data in the Ge V/TeV 
energy region during the flare favors the "low-IR" EBL model of Kneiske et al.l (120041 ). while 
a simple one-zone SSC model with a single power-law electron spectrum at injection is not 
able to reproduce the spectral shape at the highest energies (above 6 TeV). 
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Table 1. Best-Fit Parameters in the SSC Model 



Flux Level 


'Jmax 


h 


B(G) 


R (cm) 


5 


s 


Long-term 


3 x 10 6 


1.8 x 1(T 5 


0.07 


3 x 10 16 


12 


1.95 


Flare 1 


8 x 10 6 


8 x 1(T 5 


0.10 


1 x 10 16 


10 


1.6 
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Fig. 1. — Light curves of Mrk 501 with 30-day bins. The vertical bars represent the \-o 
uncertainties. The horizontal dashed lines and the legends (for all curves) show the results 
of a fit with a constant value to the data set. 
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MJD 

Fig. 2. — Light curves of Mrk 501 with 1-day bins obtained with the Swift/BAT instrument 
at 15-50 keV. The vertical bars represent the l-a uncertainties. The vertical lines show the 
start and the end of flare 1. The horizontal dashed line is the same as those in Fig.l. 




Fig. 3. — Significance map of the Mrk 501 region: (a) shows the statistical significance 
in standard deviations in the period from MJD 54683 to 55850; (b) shows the statistical 
significance of the flaring period from MJD 55851 to 55887. 
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Fig. 4. — Significance map of the Mrk 501 region during the flaring period from MJD 55851 
to 55887. Events with N pad > 500 are used, whose corresponding median energy is 8.4 TeV. 



-18 - 




Fig. 5. — Spectra of Mrk 501 during the flaring period from MJD 55851 to 55887. The 
spectrum from 0.1 GeV to 300 GeV is observed by Fermi-LAT. The solid line is the fit to 
the Fermi-LAT data using a power law. The spectrum observed by ARGO-YBJ has been 
corrected for EBL absorption using four different EBL models. Details can be found in the 
text. 
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Fig. 6. — Spectral energy distribution of Mrk 501. The solid lines show the best fits to the 
data with a one- zone SSC model, with the best-fit parameters listed in Table 1. The black 
curve corresponds to a SED model describing the long-term averaged data, while the red 
curve describes the flaring data. 



